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Abstract 

The Unitary Correlation Operator Method (UCOM) is applied to real- 
istic potentials. The effects of tensor correlations are investigated. The 
resulting phase shift equivalent correlated interactions are used in the 
no-core shell model for light nuclei and for mean-field calculations in 
the Fermionic Molecular Dynamics model for nuclei up to mass A = 48. 

1 Motivation 

In principle the interaction among hadrons should be derived from Quantum 
Chromo Dynamics (QCD), the fundamental theory of the strong interaction. 
However a description of the atomic nucleus in terms of QCD is still not 
in sight. Therefore one starts with so called realistic nucleon-nucleon (NN) 
potentials that reproduce phase shifts up to Eiab = 300MeV and the deuteron 
properties. This procedure is not unique as the NN-interaction is momentum 
dependent. Furthermore, exact calculations of the three- and four-nucleon 
system show that realistic NN-potentials are not sufficient and genuine three- 
body forces are needed to describe the binding energies. 

Our aim is to use realistic NN-potentials in nuclear structure calculations 
where the many-body Hilbert space is spanned by Slater determinants. Al- 
though those form in principle a complete basis they are very badly suited to 
describe the correlations induced by the repulsive core and the strong tensor 
part in the NN-potential V. A large scale shell model calculation with the 
bare interaction for 3 He needs already 50 hfl excitations, a Hilbert space of 
dimension 10 5 , to achieve convergence. 



2 Unitary Correlation Operator Method (UCOM) 

In order to avoid the large off-diagonal matrix elements that scatter to high 
lying shell-model states we propose [1-2] to first perform a unitary transfor- 
mation of the Hamiltonian 



H = C ] (T + V)C = T + H [2] + H m + 



(1) 



that incorporates the effects of the repulsive and tensor correlations in the 
sense of a pre-diagonalization. In Eq. (JIJ) T denotes the kinetic energy, 

the two-body, and the three-body part of the correlated Hamiltonian. 
With this unitary transformation the original eigenvalue problem in terms of 
the bare Hamiltonian H and many-body states | ^ n ) that include short ranged 

correlations can be rewritten in terms of a correlated Hamiltonian H and more 
simple states | \l/ n ) that are Slater determinants or superpositions of a limited 
number of those: 

H\V n )=E n \y n ) -> H\* n ) = C?HC\* n ) = E n \* n ) (2) 
The correlator C = CnC r consists of the unitary radial correlator 

C r = expi-iY,l(s{rij) Prij + h.a.)\ (3) 

and the unitary tensor correlator 

C n = exp{-i £ ^foj) ((£tfv)(&j m) + h.a.) ) . (4) 

When applied to a many-body state | \1/ ) the radial correlator C r shifts 
all particle pairs radially away from each other whenever they are too 
close, i.e. inside the range of the repulsive core. The strength function s(ry) 
controls the amount of the radial shift and is optimized to the potential under 
consideration. p ri j is the radial component of the relative momentum. The 

effect of the transformation \ ^f) — > C r | ^ ) is shown in the upper part of Fig. ^ 

(2) 

where the two-body density p s T is displayed as a function of the distance vector 

(ri—r-i) between two nucleons in 4 He. The p§ T on the l.h.s. is calculated 
with the shell-model state | (lsi/2) 4 ) that is just a product of 4 Gaussians. 
It has a maximum at zero distance which is in contradiction to the short 
ranged repulsion of the interaction. This inconsistency is removed by the 
action of the radial correlator C r that moves density out of the region where the 
potential is repulsive. The corresponding kinetic, potential and total energies 
are displayed in the lower part of the figure for three nuclei. The radially 
correlated kinetic energy ( C\,TC r ) increases somewhat compared to ( T ) but 
this is overcompensated by the gain of about —25 MeV per particle in the 
correlated potential energy. Nevertheless the nuclei are still unbound. 

The tensor correlations are induced by Cq where the tensor operator in 
the exponent (Eq. (pjj) creates shifts perpendicular to fy. The amount of the 



-2-10 1 2-2-10 1 2-2-10 1 2 

x rfml x ffm] x rfrn] 



4 Re 16q 40 C 


a /— 
\_ 







Figure 1: Upper part: two-body density Pst(^i~ ^2) °f 4 He for a pair of 
nucleons with isospin T = and parallel spins, S' = l,Ms = l. Arrow indicates 
spin direction and (x,y,z) = (n— r 2 ) relative distance vector. Lower part: 
corresponding kinetic, potential and total energies per particle of 4 He, 16 and 
40 Ca, without, with radial, and with radial and tensor correlations (Bonn-A). 



displacement depends on the spin directions and <jj of the particles relative 
to their distance vector. The operator pmj = p^ — p ri j, called orbital relative 

momentum, is the component of pij perpendicular to fy. The overall strength 

of the tensor correlations is controlled by $(r^-) and allows, for example, to 
map a purely I = deuteron wave function onto the exact one which includes 
an I = 2 component and thus all tensor correlations [1] . The application of the 
tensor correlator Cq leads to the two-body density depicted in the right hand 
contour plot of Fig. ^ One may visualize the action of Cq as a displacement 
of probability density from the 'equator' to both 'poles', where the spin of 
the 5 = 1 component of the nucleon pair defines the 'south-north' direction. 
Again this costs kinetic energy but now the many-body state is in accord with 
the tensor interaction and one gains the binding needed to end up with about 
-8 MeV per particle, see lower part of Fig. ^ 



3 Choice of correlation functions 



Like the interaction the correlators are decomposed into the four spin isospin 
channels S— 0,1; T = 0,1 and the corresponding strength functions s(r) and 
i9(r) are adjusted separately. As the repulsion is of short range the optimal 
radial shift functions can be found by minimizing the energy with respect 
to s(r) in the corresponding channel. The result depends only weakly on 
the choice of the system, a constant trial state in the two-body system gives 
practically the same s(r) as minimizing the energy of the deuteron or 4 He. 

Unlike the repulsive core the tensor part of the interaction, that is due 
to pion exchange, is of long range. Therefore, adjusting $(r) by mapping a 
pure / = trial state for the deuteron to the exact eigenstate leads to the very 
long ranged tensor correlator strength $ (r) shown in Fig. When used in 
systems with more than two particles this correlator is not useful as it induces 
large three- and more-body parts in the correlated interaction (H^\ . . .) that 
are very complicated to handle. Therefore we restrict the range of the tensor 
correlations and consider in the following the three $(r) labeled by a, (3 and 7 
shown in Fig. |21 With those we take care of the short range part of the tensor 
correlations. The long range part must then be described by the many-body 
state I ^/). The presence of the other particles will destroy at least partially 
the alignment of spins of a particle pair at larger distances (spin frustration) 
so that I \l/ ) can be again a rather simple many-body state. 
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Figure 2: Strength of tensor correlation as function of particle distance. 



4 Correlated interaction in momentum space 

In Fig. El we display in momentum representation the matrix elements of the 
bare Argonne V8' potential with those of the corresponding correlated inter- 
action H^. The left column is for the 1 5'o channel where due to S = only 
the radial correlator acts. 

It is obvious that the goal of pre-diagonalization is achieved. Beyond mo- 
mentum transfers of about 2 fm -1 the off-diagonal matrix elements calculated 
with correlated states are close to zero. Our result is in good agreement with 
Vi ow -k obtained with renormalization group methods [1,4]. In the right column 
the matrix elements between the / = and I = 2 triplet channels are shown. 
Here only the tensor components of the interaction contribute. For the corre- 
lated state we use the correlation strength # a (r) shown in Fig. El Despite the 
restricted range the correlator achieves a substantial reduction of the matrix 
elements. With the long ranged correlator $ d (r) derived from the deuteron the 
off-diagonal matrix elements vanish completely. 
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Figure 3: Matrix elements in momentum representation for the uncorrelated 
and correlated Argonne V8' potential. Left: 1 Sq channel, plane denotes zero, 
vertical range from -2 to 2 fm . Right: tensor mixing between 3 5i and 3 D 1 
channel, vertical range from -1.5 to fm~ . 



5 No-core shell model calculations 



In order to test the UCOM we use the correlated Hamiltonian in two-body 
approximation with the no-core shell model code of P. Navratil [6] and 
compare with benchmark calculations for 4 He [5]. 

In Fig. |3] we show the ground state energy as a function of the oscillator 
parameter for 0, 4, 10 and 18 excitations calculated with tensor correlators 
of range j3 and a. The correlated Hamiltonian obtained with the medium 
ranged correlator (3 gives with a single Slater determinant (OhQ) a minimum in 
energy at almost the exact energy (horizontal line). Admixing more and more 
configurations lowers the energy and convergence is reached at about 8hQ. 
One should keep in mind that the bare interaction needs many-body spaces 
that are orders of magnitude larger. For example at hfl = 25 MeV the binding 
energy calculated with excitations up to 18 hQ included is only —6 MeV. If we 
use the short-ranged tensor correlator a the optimal single Slater determinant 
is 7 MeV above the exact eigenvalue. As can be seen in Fig. H] convergence is 
not as fast but the overbinding is reduced to about 2 MeV. This difference to 
the reference calculations is attributed to the missing contributions from 

and H®. 
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Figure 4: Ground state energy of 4 He calculated with Argonne V8' as function 
of harmonic oscillator parameter hfl and number of major shells included. 
L.h.s.: Hamiltonian correlated with tensor correlator of range /3, r.h.s.: shorter 
range a. Horizontal line: benchmark calculations [5]. 



6 Mean-field calculations in FMD basis 



For heavier nuclei no-core shell model calculations are not feasible and we 
perform mean-field calculations in the framework of the Fermionic Molecular 
Dynamics (FMD) model [3]. The trial state in this approach consists of a 
single Slater determinant with single-particle states that are parametrized as 
a single or a sum of two Gaussians. As Hamiltonian we use the correlated 
Bonn-A interaction that is complemented by a momentum-dependent two- 
body correction term. This correction term has to simulate the effect of the 
missing three- and more-body terms of the correlated interaction and the effect 
of the genuine three-body forces. It is adjusted to reproduce the binding 
energies and radii of the doubly-magic nuclei 4 He, 16 and 40 Ca. Minimizing 
the binding energy with respect to the parameters of the single-particle states 
yields intrinsic states. The resulting one-body densities are shown for some 
nuclei in Fig. E] 

Besides the spherical nuclei 4 He, 16 and 40 Ca we obtain also an axial 
20 Ne and a triaxial 24 Mg. The energies of the intrinsic states are shown in 
Fig. |U1 For the p-shell we achieve a substantial improvement by using single- 




Figure 5: One-body densities of 4 He, 16 0, 40 Ca and 12 C, 20 Ne, 24 Mg. 
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Figure 6: Deviation of intrinsic energies from experimental ground state ener- 
gies. Insert: two Gaussians per single particle state. 

particle states with two Gaussians per nucleon. For the doubly-magic nuclei 
the spherical intrinsic states can be regarded as a good approximation of the 
ground states. On the other hand nuclei between closed shells are intrinsically 
deformed so that a projection on good angular momentum should be performed 
before comparing with experimental binding energies. 

Configuration mixing calculations taking into account rotations and vibra- 
tions are under way and should allow a description of the low lying spectra. 
First results are promising. 

We are grateful to Petr Navratil for providing his no-core shell model code 
used in the 4 He calculations. 
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